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Abstract 

In Liang et al (2009), the current authors demonstrated that BSDEs 
can be reformulated as functional differential equations, and as an ap- 
plication, they solved BSDEs on general filtered probability spaces. In 
this paper the authors continue the study of functional differential equa- 
tions and demonstrate how such approach can be used to solve FBSDEs. 
By this approach the equations can be solved in one direction altogether 
rather than in a forward and backward way. The solutions of FBSDEs are 
then employed to construct the weak solutions to a class of BSDE systems 
(not necessarily scalar) with quadratic growth, by a nonlinear version of 
Girsanov's transformation. As the solving procedure is constructive, the 
authors not only obtain the existence and uniqueness theorem, but also re- 
ally work out the solutions to such class of BSDE systems with quadratic 
growth. Finally an optimal portfolio problem in incomplete markets is 
solved based on the functional differential equation approach and the non- 
linear Girsanov's transformation. 
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1 Introduction 



Backward stochastic differential equations (BSDE) provide a new perspec- 
tive to look at the infinitesimal behavior of Markov processes, and they have 
been found intrinsically linked to a class of nonlinear partial differential equa- 
tions (PDEs). In a fundamental paper by Pardoux sand Peng [33], they solved 
a class of nonlinear BSDEs with Lipschitz drivers for the first time. An intrinsic 
connection to nonlinear PDEs, which is now well known as nonlinear Feynman- 
Kac formula, was later established by Peng [37J and Pardoux and Peng [35] . Its 
applications in finance were discovered by Duffle and Epstein |13| and El Karoui 
et al [IS]. For a more comprehensive review of the BSDE theory, we refer to 
[T4] . [T5] and [41] and reference therein. 

In Liang et al [26], the current authors demonstrate that BSDEs can be 
reformulated as functional differential equations. As an application we can 
solve BSDEs on general filtered probability spaces, and in particular without 
the requirement of martingale representation. In this paper we try to apply such 
functional differential equation approach, or functional approach for short, to 
forward backward stochastic differential equations (FBSDEs), and demonstrate 
how it can be used to solve a financial optimal portfolio problem in incomplete 
markets. 

Let us recall such idea, which is the first main ingredient of the paper. 
Let (f2, J 7 , Tu P) be a filtered probability space satisfying the usual conditions. 
Given a semimartingale [o,t] with the following decomposition: 

Y t =M t -V t , fortG[0,T], 

where M is local martingale and V is a finite variation process, if we further 
know the terminal data of Y, say Yr = £ for some J-r-measurable random 
variable £, then we also have £ = My — Vt- Therefore 

f Y t = E[Z + V T \Ft]-V t , 

\ M t = E[£ + V T \Tt], for t G [0, T]. 1 ' ' 

Before we apply the above relationship to a specific FBSDE setting, 

let us first look at it from potential theory point of view. For a given domain 
D C R d , a real- valued function / is called superharmonic in D if 

/ fix + y)o~ r (dy) < f(x), for x G D and r < dist{x, dD), 

JdB(0,r) 

where oy is the measure on the surface of the ball dB(0, r), normalized to have 
the total mass 1. / is called harmonic if furthermore the equality holds. If / 
is superharmonic in D, then there exists a unique positive Borel measure on D 
such that the following Riesz decomposition holds: 

f(x) = h(x) + Gd/i(z), for x G D, 

where h is harmonic in D and Got 1 is the Green's potential of fi on D, i.e. 
Gol^ix) = J D Gix, y)fi(dy) with G(x,y) as the Green's function of the Laplace 
equation in R d . Given the boundary data of such superharmonic function /, 
by the above Riesz decomposition, the Green's potential GdH is often used to 
study /. 
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The probabilistic counterpart of the above Riesz decomposition is the Doob- 
Meyer decomposition. For a supermartingale (Yt)te[o,T] with Cadlag sample 
paths, the Doob-Meyer decomposition says that there exists a unique increasing 
predictably measurable process V starting from Vq = such that M defined by 
Mt =Y t + Vt for t G [0, T] is a martingale. The above relationship we just 
established tells us exactly the same thing as in the potential theory: given the 
terminal data of a supermartingale (which is a scmimartingale) defined on [0, T], 
we can study such supermartingale by investigating the increasing predictably 
measurable process V. 

Now we apply the above relationship (|l.ip to a specific FBSDE, which will be 
served as an auxiliary equation later. For any given filtered probability space 
(fi, J 7 , Tt , Q) satisfying the usual conditions, with a d-dimensional Brownian 
motion B = (B , • • • , B d )* , (The superscript * denotes the matrix transposi- 
tion) 



dX t = f(t,Y t ,Z t )dt + dB t , 
I =i€ R d , 

dY t = -h(t,Y t ,Z t )dt + Z t dB u 
{ Y T = (t>(X T ). 



(1.2) 



with the coefficients / : [0, T] x R n x R nxd -> R d , h : [0, T] x R n x R nxd -> i?" 
and : R d —> R n satisfying Condition [T] to be introduced later. 

To solve FBSDE (|1.2[) , for r G [0, T], we consider the following functional 
differential equation on [r, T] : 



V t = J h(s,Y{V,X) s ,Z(V 1 X) s )ds 
together with the forward process X: 

X t =x+ f f(s,Y(V,X) a ,Z(V,X) s )ds+ I dE 



where 



f Y(V, X) t = E[(j>(X T ) + V T \T t } - V t , 

T 

Z(V,X) s dB s = 4>(X T ) +V T - E[4>{X t ) + Vt\F t ]- 



(1.3) 



(1.4) 



(1.5) 



If we can solve (V, X) for functional differential equation (11.311 (jl.4H with r = 
0, then by Lcmma[5]in Section 2.2, (Y(V,X), Z(V,X),X) will provide us with 
the solution to (|1.2p . For the notation's simplicity, let us denote = (V,X)* 
and F = (h, /)*. Then ([OJ (JOJ are simplified to: 



% = Xx + J F(s,Y(^) 9 ,Z(^) s )ds + J X idB s 



(1.6) 



where G R n+d with the first n components being and next d components 
being x G R d . Note that by functional differential equation (|1.6[) . we can solve all 
the components (Y, Z, X) in one direction altogether. We consider the solutions 
to (|1.6[) in the following space: 
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• C([0, T}; R n+d ): the space of continuous and J^-adapted processes {^t)te[o.T] 
valued in R n+d such that sup tg r r i \^t\ £ £ 2 (f2, J-"t, Q) and endowed with 
the following norm: 

ll*Hc[o,r] = Je sup |* t | 2 . 
V *e[o,T] 

In Section 2.2 we will mainly solve functional differential equation (|1.6p and 
prove the following theorem: 

Theorem 1 // the coefficients satisfy Condition [7] then there exists a unique 
solution ^ £ C([0, T}; R n+d ) to functional differential equation U.6}) . 

The second main ingredient of this paper is a nonlinear version of Girsanov's 
transformation, which is employed to connect the above FBSDE (|1.2j) and a class 
of BSDE systems with quadratic growth. Namely wc will use strong solutions of 
FBSDEs to construct weak solutions to a class of BSDE systems with quadratic 
growth. 

In the PDE theory, if the nonlinear terms in equations have at most quadratic 
growth with respect to the gradient of solutions, the nature of the equations com- 
pletely change. In the BSDE theory, there is a class of BSDEs with quadratic 
growth corresponding to such PDEs, and they are usually called quadratic BS- 
DEs. The study of quadratic BSDEs was initialized by Kobylanski [33J using 
the idea of the Cole-Hopf transformation adapted from the PDE theory. Her 
result was substantially developed and generalized by Briand and Hu [5] and 
[5], where they extended to the equations with the unbounded terminal data 
and with the convex driver. On the other hand, Quadratic BSDEs have found 
a lot of applications in finance. For example, they appear naturally when one 
wants to derive the value function for the maximization of expected utility, use 
indifference pricing idea to hedge a contingent claim written on nontradeable 
underlying assets, or consider the risk measure. See for example [1] [T7j [H] [35] 
[33] and [39]. 

In this paper we mainly consider the following quadratic BSDE system (not 
necessarily scalar): 

dY t = -h(t, Y t ,Z t )dt - Z t f(t, Y t ,Z t )dt + Z t dWt, 

(1.7) 

Y T = <j)(W T ) 

where W = (W 1 , ■ ■ ■ , W d )* is a d-dimcnsional Brownian mtion starting from 
x £ R d . The coefficients h, f and <f> arc supposed to satisfy the following 
condition: 

Condition 1 All the coefficients h : [0,T] x R n x R nxd -> R n , f : [0, T] x 

R n x R nxd — > R d and <j> : R d — > R n are continuous. Moreover h, f and <p are 
Lipschitz continuous, i.e. 

\h(t, y, z) - h(t, y, z)\ < d{\y -y\ + \z- z\), 
\f(t,y, z )-f{t,y,z)\<C x {\y-y\ + \z-z\), 
\4>{x)-4>{x)\<C 2 \x-x\, 
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and <fi is uniformly bounded, 

sup \<j>(x)\ < M, 
xeR d 

forte [0,T], y,y £ R n , z,z £ R nxd and x,x£ R d . 

Because of the terms with the coefficient /, the equations have at most 
quadratic growth, i.e. there exists a constant C3 such that for any y £ R n and 

Z £ R nxd , 

\zf(t,y,z)\<C 3 \z\(t+\y\ + \z\). 

The quadratic growth term in (|f ,7|) is more special than the usual one consid- 
ered in the literature. However this special structure is enough to cover the 
most examples of quadratic BSDEs known in finance, at least with some extra 
conditions added. We will consider one specific example from optimal portfolio 
problems in Section 4. Moreover, all the excising results of quadratic BSDEs 
are only for the case n = 1. The current paper seems to be the first attempt to 
consider the quadratic BSDE systems. 

On the other hand, one may wonder why the terminal data has the special 
form <P(Wt)- This is only for the presentation's simplicity. The whole paper's 
results can be extended without difficulty to the case 4>(Xt) where X is driven 
by stochastic differential equations (SDEs): 

dX\ = X\{b\dt + a l t dW t ), for i = 1, • • • , m. 

with the coefficients b z and a 1 satisfying certain regularity conditions. 

To solve (|f ,7p . we will pursue another direction different from the existing 
method for quadratic BSDEs. Namely we don't use the Cole-Hopf transforma- 
tion at all and don't assume the underlying probability space and Brownian 
motion as any given; instead we consider weak solutions of quadratic BSDEs. 

Before presenting the definition of weak solutions to (|f .7[) , let us mention 
some already existing work about weak solutions. One of the first attempts to 
introduce the weak solutions for BSDEs was Buckdahn et al [7J, and Buckdahn 
and Engelbert [5] further proved the uniqueness of their weak solutions. However 
the driver of their BSDE does not evolve the martingale representation part Z . 
On the other hand, the notion of weak solutions for FBSDEs was introduced by 
Antonelli and Ma [5] and further developed by Delarue and Guatteri [TJ], and 
by Ma et al [30] and Ma and Zhang [23] who employed the martingale problem 
approach. 

Definition 2 A weak solution to BSDE |1.7| j is a triple (Vl, F ,~P X ) , {Tt} and 
{Y,Z FX ,W) such that 

(1) (Q, J 7 , ~P X ) is a complete probability space with the filtration {Ft} satis- 
fying the usual conditions; 

(2) under such filtered probability space, Y, Z F and W are Tt-adapted, and 
Y is a special semimartingale, Z p is the density representation ofY under P x , 
and W is a Brownian motion starting from P x (Wq = x) = 1; 

(3) The increments {W u — Wt ■ t < u < T} must be independent of the 
a-algebra J-" t ; 

(4) the following integral equation satisfies: 

Y t = q>{W T ) + J h(s,Y s ,zf)ds + J Zff(s,Y s ,Zf)ds- J zfd\Y s . 

(1.8) 
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Remark 3 In this paper, the density representation Z means it is the density 
representation for the martingale part of the special semimartingale Y , and we 
use the superscript P x to emphasize the dependency of the density representation 
on the probability measure P 1 . 

Remark 4 Our definition of weak solutions is more related to Buckdahn et 
al The filtration {J 7 *} plays an important role here. If Ft = > *- e - the 
filtration is generated by the Brownian motion W augmented by the P x -null sets 
in F, the solution turns to be a strong solution. In Ma and Zhang \31}/ . such 
solution is also called a semi-strong solution. Actually the smallest filtration for 
weak solutions is the filtration {F^ ,Y,Z } generated by W, Y, Z and satisfying the 
usual conditions. 

Remark 5 Condition (3) automatically holds given the Brownian motion W 
with the filtration {Ft}- In fact such condition simply means {Ft} consists, 
additionally to {F^}, only of independent experiments. In Buckdahn et al 
such condition is formulated in terms of martingales, i.e. any F^ -martingale 
must be an Ft-martingale. In Kurtz 124-1. such kind of condition is called the 
compatibility constraint. (3) is extremely useful when we want to identify weak 
solutions are strong solutions. 

For the notation's simplicity, we will suppress the superscript x of P x from 
now on if no confusion may arise. Now we describe our idea formally before 
presenting the existence and uniqueness theorem of BSDE p.7[) . The basic 
idea is using the strong solution of FBSDE to construct the weak solution to 
quadratic BSDE. Let us start with a Brownian motion family B on (f2, F, Q) 
with the filtration {Ft} satisfying the usual conditions, and consider FBSDE 

Suppose FBSDE (|1.2p admits a unique solution (X,Y,Z Q ). Then we define 
a new probability measure P by 

£-™ 

where S (N) is the Doleans-Dade exponential of N with 

N = — [ (f(s,Y s ,Z?),dB s ) d 
Jo 

where (•, -}d denotes the inner product in R d . Under the new probability measure 
P, by the Girsanov's theorem, B has the following decomposition: 

B= (B-[B,N]) + [B,N], 

= (b + J f(s,Y 3 ,Z?)ds) - J f{s,Y s ,Zf)ds 

where B — [B, N] = B + J Q f(s, Y s , Z^)ds is a martingale under P, and further- 
more by the Levy's characterization, it is in fact a Brownian motion under P. 
We further define W by W = x + B - [B, N]. 

Under the probability measure P and with the Brownian motion W, let's 
rewrite the backward equation in FBSDE (|1.2p : 

dY t - -h(t, Y u Z?)dt - Z?f(t, Y u Z?)dt + Z?dW t 



G 



with Y T = (j)(W T ). If we can prove Z Q = Z p , then triple (fi,.F,P), {F t } and 
(Y, Z p , W) is just one weak solution we want to find. 

There are mainly three steps needed to be verified for the above solving 
procedure. The first step is about the invariant property of the density repre- 
sentation under the change of probability measure, i.e. Z® = Z p ; The second 
step is of course the solvability of FBSDE (|1.2|) ; the last step is the Doleans- 
Dade exponential S{N) must be a uniform- integrable martingale to guarantee 
P is a probability measure. As long as the above three steps are verified, we 
have the following theorem: 

Theorem 6 If the coefficients satisfy Conditional^ then there exists at least one 
weak solution (fi,.F,P) 3 {T t } and (Y,Z P 1 W) to BSDE [Tty . 

The paper is organized as follows: In section 2 we verify the above three steps 
and prove Theorem [SJ while in Section 3 the uniqueness and the connection 
between weak solutions and strong solutions are discussed. Finally we apply 
our method to an optimal portfolio problem in incomplete markets in Section 
4. 

2 Weak solutions and existence 

2.1 Invariant property of density representation 

The following lemma is almost trivial but crucial to our results, which states 
that the density representation of a special semimartingale is invariant under 
the equivalent change of probability measure. This observation is firstly made 
in Liang et al [27], and we recall it here for completeness. 

Lemma 7 Let B be a Brownian motion on (Q, J 7 , Q) with the filtration {J^t} 
satisfying the usual conditions. Let Z^ be the density representation of a special 
semimartingale Y under Q. If define an equivalent probability measure P by 
= S'(N) for some uniform integrable martingale S'(N), then Zf = for 
a.e. t £ [0, T], a.s. 

Proof. Under the probability measure Q, Y has the canonical decomposi- 
tion Y = M — V with M being a local martingale and V being a finite variation 
process, and moreover, M admits the martingale representation: 

M t -M = f ZfdB s , for t e [0,T], a.s. (2.1) 
Jo 

for some predictably measurable process Z®. 

By Girsanov's theorem, Y is still a special semimartingale under P but with 
the canonical decomposition Y = M—V, where M = M—[M, N] is a martingale, 
and V = M — Y is a finite variation process. We also have B = B — [B, N] as 
a Brownian motion under P. Hence under P, (|2.1[) becomes 

M t - M + [M, N] t = f ZfdB s + f Z?d[B, N] s , for t e [0, T], a.s.. 
Jo Jo 
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By identifying the martingale parts and finite variation parts of the above equal- 
ity, we must have 



M t -M 




o 



But on the other hand under P we also have 




ZfdB s , for t G [0,T], a.s. 



Jo 

for some predictably measurable process Z F , 



so 




f-Z?\ 2 ds = 0, a.s., 



which proves the claim. ■ 

Since the density representation usually determines the hedging (or repli- 
cating) strategy in finance, a direct consequence of Lemma [7J is that hedging 
strategy is independent of the choices of the equivalent (martingale) probability 
measures. Due to Lemma [7J we will not emphasize the dependency of the den- 
sity representation on the probability measure, and simply write it as Z from 
now on. 

2.2 Functional approach to FBSDEs 

The study of FBSDEs was initiated by Antonelli [1], and this subject was 
further developed in J22] [H] [3H] [3H] [10] and especially the monograph [55] by 
Ma and Yong. However most of them either solved the equations locally or as- 
sumed some regularity on the coefficients (e.g. smoothness and monotonicity). 
Recently Delarue [TU] solved FBSDEs globally with Lipschitz continuous as- 
sumptions on the coefficients by combining the method of contraction mapping 
and the four-step scheme of FBSDEs. 

In this subsection we try to reformulate FBSDE (|1.2j) as functional differen- 
tial equation (|1.6|) and solve such functional differential equation instead. The 
approach may benefit especially numerical solutions of FBSDEs, because a usual 
obstacle to numerically solve FBSDEs is one needs to solve (Y, Z) backwards and 
X forwards at the same time. By introducing a functional differential equation, 
we can solve all the components (Y, Z, X) in one direction altogether. 

We first establish the equivalence between FBSDE (|1.2[) and functional dif- 
ferential equation (|1.6p . Besides the space C([0, Tj; R n ) we further introduce the 
following space: 

• H 2 ([0,T}; R n ): the space of predictably measurable processes endowed 
with the norm: 



Lemma 8 FBSDE fO)) admits a unique solution (Y,Z,X) e C([0,T];R n ) x 
-ff 2 ([0, T]\ R nxd ) x C([0,T]; R d ) if and only if functional differential equation 
il.6]) admits a unique solution \I/ e C([0, T}; R n+d ), and therefore by Theorem 
1, FBSDE admits a unique solution. 
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Proof. Suppose (Y, Z, X) is the unique solution to FBSDE (TO]). Since Y is 
a special (continuous) semimartingale, it admits the following canonical decom- 
position: Y t = Mt — Vt for t <G [0, T], where M is a continuous local martingale 
and V is a continuous finite variation process with Vo = 0. Furthermore by the 
martingale representation M t = Mq + j Q Z s dB s for t € [0, T], we have 



F 4 = M + / Z s dB s - V t , for t e [0, T] , 
Jo 



(2.2) 



from which we obtain the relationship (jl.5D . The backward equation in FBSDE 
(11.21) becomes 



Mo + / Z s dB s -V t = <t>{X T ) + / f(s, Y s ,Z s )ds 



Z x dB s 



By taking conditional expectation with respect to Tt on both sides: 



M 



Z s dB s - V t 



n 



E[<t>(X T )\Ft] + E 
E[<t>{X T )\F t ]+E 



J f{s,Ys,Zs)ds\Ft 

f f{s,Y s ,Z s )ds\F t - [ f(s,Y s ,Z s )ds. 
Jo Jo 



By the uniqueness of the canonical decomposition of Y, and by identifying the 
martingale part and finite variation part of both sides, V must satisfy (II. 3|) . 

On the other hand, if \1/ is the unique solution to functional differential 
equation (11.61) . the relationship (|1.5[) can be rewritten as 

Y(V,X) t = E[<j>{X T ) + V T \T t ] - V t 

= E[<j>(X T )+V T \Ft}- f f(s,Y(V,X)s,Z(V,X)s)ds, 
Jo 

together with 

J Z s dBs=^X T ) + V T -E[^X T ) + V T \J : t], 

from which we deduce {Y(V, X), Z(V, X), X) must satisfy ([P]) . ■ 

The rest of this subsection is devoted to the proof of Theorem 1. 
Lemma 9 If the coefficients satisfy Condition [IJ and r satisfies 

then functional differential equation il.b)) admits a unique solution £ C([t, T]; R n+d ). 
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Proof. The mapping defined by (|1.6[) is denoted by L. We will first show 
that L : C([t, T}; R n+d ) -> C([t, T]; In fact for * e C([t, T]; 

im*)lle[T,T] 



|F(s,y(*) s ,Z(*) a )|ds 



'E sup | / xi^ s | 2 

T<t<T Jt 



< \x\ + Vt^ 



E 



£ ' \F(8,Y(9) a ,Z(*).)\*<kj + 2^E\ j\idB a \* 



< \x\ + (Ci\/T^7+2rf)W / (s 2 Vl)d 



CxVt^JeJ^ \Y(v) B \*ds + dy/r^Fd e J |z(*) s | 2 d s . 



Note that 



'E J \Y(^) s \ 2 ds 



<\j E J r {E[cf>(X T )\T s }} 2 d S + J E {E{V T \F s ]}2ds + JE j \V s \*ds 



< 



E\cj>{X T )\Hs 



E\V T \ 2 ds- 



E\V s \ 2 ds 



< v^7(2 + C 2 )||*|| c[T)T] , 
and by Ito's isometry, 



'eJ |Z(*) s | 2 ds 



Z(V) s dB s 



< s/E\ct>{X T )\ 2 + v/£|VH* + JE{E[<f>(X T )\F T ]} 



1 E {E[V T \F T \} 2 
< (2 + 2C 2 )||*|| c[T , T] . 

Therefore | |L(*)| | c[t>t] < oo. Similarly for e C([r, T]; R n+d ), we have 

||£(*)-£(*')|| C[T , T] 

< CiVr^(v / r^(2 + c 2 ) + 2 + 2C 2 ) ||*-*'|| c[r , T ] 

< 2ll* _ *'llc[T,n 



by the condition on r. Hence L defined by (|1.6[) is a contraction mapping on 
C([r,T];R n + d ). m 
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Based on Lemma [§J we next extend to the global solution on [0,T]. To do 
this we pursue the bounded solutions for FBSDE (|1.2p . First by the Markov 
property, there exists a Borel-measurable function $ such that Y t = <$>(t,X t ). 
By checking the proof for Lemma |H1 the crucial step to extend to the global 
solution of (|1.6[) is that one needs a uniform estimate for the gradient of <1>. We 
recall a regularity result from Delarue . 

Lemma 10 (Delarue Ulf ) Under Condition^ on the coefficients, there exists a 
Borel measurable $ such that Y t = X t ). Moreover there exists a constant C4 
depending on the Lipschitz constants C\ and C2, the bound M of the terminal 
data, the dimension n and d, and the terminal time T such that 

|$(t,z)|, |V x $(t,i)| < Ci, for (t,x) S [0,T] x R d . 

Based on such constant C4, we make a partition of [0,T] by tt : — to < 
t± < • ■ ■ < tjsr = T with the mesh = maxi<i<jy \U — such that 

1 Al. 



8C 1 (i + a 4 ) 



We start with \t N ^t N \ and consider *(iV) = (V(N), X(N) tN -^ x )* such 
that 



with 



*(N) t = Xx + F(s,Y(N) s ,Z(N) s )ds+ xidB. 



Y(N) t = E[^{X(N) t ^- 1 ' x ) + V{N) T \F t ] - V(N) t , 



Z{N) s dB s = 4>(X(N)^- UX ) + V(N) T 



where we used the superscripts (4/v-i, x) to emphasize X(N) tN - 1,x starting from 
X(N) t ™2l ,X — x. By Lemma IH1 there exists a unique solution: 

*(JV) = (V(N),X(N) tN -^ x )* eC([t N -i,T];R n+d ) 

and we also get (Y(N), Z(N)). Moreover, there exists a Borel-measurable func- 
tion $jv-i such that Y(N) tN _ 1 = $jv-i(^w-i, x) and by Lemma [Till 

|V x $iv-i(tjv-i,a;)| < C 4 forarei^. 

In general on [U-i,U] for 1 < i < N- 1, consider *(z) = ^(i)' 1 - 1 ^)* 
such that 

*(*)* = Xx+ / F(s,Y(i) s ,Z(i) s )ds+ [ X idB s 



with 



r(i) t = + ^(OtJ^] - v(i) t 

Z{i) s dB s = Mt t ,X(i)%- uX ) + V{i) u 
ti-1 

•\ti— 1 
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By Lemma [5] again, there exists a unique solution: 

and we get (Y(i), -Z(i)) as well. Moreover there exists a Borel-measurable func- 
tion 3>i_i such that Y(i) ti _ 1 = 3>j_i(£s_i, a;) an d by Lemma ITUl 

|V x $i-i(ti_i,a;)| < C* 4 forxei? d . 

Of course (V(i), X(i) ti_1,a: ) for 1 < i < JV are not the real solutions to (|1.6[) 
on the corresponding time interval because they start from 

We need to shift the paths of (V(i),X(i) ti - 1 > x ) accordingly in order to match 
the starting points for the solutions to (|1.6[) on each time interval [U-i, ti]. 

Lemma 11 If the coefficients satisfy Condition^ then the global solution ^ = 
(V, X)* to U.6]) is constructed as follows: for 1 < i < JV, 



Vt = V(i) t + V(j) tj for U-i <t<U 



where we follow the convention Y^°j=i = and 



X, 



X{2) 



for t Q <t<tx] 
fort! <t<t 2 ; 



{ X(N)t 1 fort N -!<t<t N 

with (Y, Z) being constructed as (Y t ,Z t ) = (Y(i) t , Z(i) t ) for <t< ti. 

Proof. We only need to show \I> = (V, X)* with (Y, Z) satisfying (|1.6[) for r = 0. 
In fact for i S [tjv-i,tjv], by the definition of (V, X), 



JV-l 



JV-l 



V* - = v(jv) t + Y, v(jh - - E y w* 



/i(s,y(jv) s ,z(jv) s )(j s , 



and 



so 



X t - X tN _, = X(N) 



tw-i-.Xt 



f{s,Y(N) a ,Z(N) B )ds 



dB„ 



tjv-i 



/ F(s,Y(N) s> Z(N s ))ds+ f X idB s , 
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where 



Y(N) t =E[<P (x(N)*f - 1 ' + ^(JV)r|Ji] - V(W)t 

N-l 



= E 



3=1 



and 



JV-l 
3=1 

.B[0(x T ) + i/ T |j- t ]-y t , 



Z(N) s dB s = ct>(X(N) t r ?- 1 ' XtN - 1 ) + V{N) T 



E 



0(X(N) 



tN — 1 >Xt 



+ V{N) T \E ti 



N-l 



4>{x T ) + v(n) t + v ^)t 



3 = 1 



-E 



N-l 



4>(X T ) + V{N) T + Y,V{j) tj \F tN _ x 

3 = 1 

= 4>{X T ) + V T - E[4>(X T ) + VrlFts-r]- 

Hence (V,X) with {Y, Z) defined in the lemma satisfy (|1.6I) on [tjv-i,tjv]- 

In general for 1 < i < N — 1, by the backward induction, it is easy to verify 
* t = (y t ,X { )* with (y t , Z t ) also satisfy (HBJ for t £ [tj-i,**]. ■ 

2.3 Uniform integrability of stochastic exponential 

In this subsection we will verify the Doleans-Dade exponential S'(N) is a 
uniform-integrable martingale. To prove this we need an appropriate martingale 
space. It turns out the martingale space we need is bounded-mean-oscillation 
(y37\/0)-martingale space. As the BMO-martingale theory is already quite 
standard in the quadratic BSDE study (e.g. [3] and [H]), we only recall some 
basic facts that are necessary in the following. For the further details and proofs, 
we refer to He et al [IB] , 

Let M be a continuous local martingale on [0,T]. For p > 1, define the 
martingale space W equipped with the norm ||M||?#> = E{[M, M}^/ 2 } 1 ^. For 
p > 1, W is the dual space of H q with q being the conjugate of p, i.e. 1/p+l/q = 
1. However for p = 1, the dual space of H 1 is strictly larger than the 
class of all martingales with bounded quadratic variation. By the Fcfferman's 
inequality, the dual of H 1 is in fact BMO2, which is the subspace of T-L 2 and 
such that there exists a constant C5, 

E {\M T - M T \ 2 \E T } <C 2 

for any stopping time r < T, and C5 is defined to be the BA/CVnorm. Similarly 
we can also define i?MO p -space for any p > 1, which are equivalent to each 
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other. In fact for p > 1 and M £ T-L 2 , there exists a constant Ce(p) depending 
on p such that 



\\M\\ BM o 1 < \\M\\ B mo < C 6 (p)\\M\ 



BA/Oi 



by the Jensen's inequality and the John-Nirenberg inequality respectively So 
from now on we will simply write BMO without specifying p. 

Lemma 12 If the coefficients satisfy Condition[l\ then 

N = - [ (f(s,Y 3 ,Z s ),dB s ) d 



is a BMO -martingale under Q, and therefore the Doleans-Dade exponential 
S'(N) is a uniform-integrable martingale under Q. 

Proof. For any stopping time r < T, by Ito's isometry and the linear growth 
condition \ f(t, y,z)\ < C\(t + \y\ + \z\), we obtain 

sup^[|^ T -iV T | 2 |^ r ] 



sup E 



J \f(t,Y tl Z t )\ 2 dt\F T 













< CfT 3 + 3C 1 2 sup£ 


J \Y t \ 2 dt\F T 


+ 3Cfsu-pE 


J \Z t \ 2 dt\F T 


, a.s 


T 




T 







(2.3) 

Since Y is uniformly bounded, we only need to control the last term of (|2.3D . By 
applying Ito's formula to (Y t ) 2 from r to T and taking conditional expectation 
on J- T , wc obtain 



\YrV+E 



\z t \ 2 dt \ j; 



E[(/)(XT) 2 \J r T }+2E 



(Y t ,h(t,Y t ,Z t )) n dt\T T 



< M 2 + X 2 E 

< M 2 + \ 2 E 

3C 2 



A 2 



■E 



\Y t \ 2 dt\F T 



\Y t \ 2 dt\T T 



\Z t \ 2 dt\F T 



\h(t,Y u Z t )\ 2 dt\T T 



C 2 T 3 iC 2 



E 



\YA 2 dt\F T 



a.s., 



where we used the elementary inequality 2ab < X 2 + b 2 /X 2 . By choosing A large 
enough such that 1 — 3C 2 /A 2 > 0, and by the uniform boundedness of Y , there 
exists a constant CV such that 



axvpE 



\Z t \ 2 dt\T T 



< Cy, a.s. 



and the conclusion follows by plugging the above estimate into 
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3 Uniqueness and strong solutions for n = 1 



As in the classical SDE theory, there are also several notions of uniqueness 
for BSDEs as well. In this section we discuss the uniqueness of BSDE (|1.7p . Wc 
assume the following condition on the coefficients: 

Condition 2 Condition]]] is assumed to be satisfied. Moreover n = 1, i.e. 
BSDE [Hfy is a scalar BSDE; F = F(t,z) with F = (h,f)*, i.e. both of the 
coefficients h and f only depend on t and z; and f J = f 3 (t, z J ) for j = 1, • • • , d, 
i.e. there is no mixture terms of z in f . 

The weak solution can be regarded as probability distribution on the sample 
path space, so we will specify the sample path space of (|1.7[) firstly. By W m we 
denote the space of continuous functions C([0, T]; R m ). Define the coordinate 
mapping X t : W m -> R m by 

X t (x)=x t , forxeW™, 

and on W m , define the following cr-algebras: 

B x = a(x s :s<t); Bf = a(x u - x t : t < u < T); 

and B x = V te[0iT] S t x . Obviously we have the relationship B x = Bf V Bf 
for any t £ [0, T]. By Definition [2J Y and W must be continuous. However it 
is not obvious at all that the density representation Z has any path regularity. 
Fortunately under Condition [2 Imkeller and Dos Reis [25] already did this work 
for us, which states that there is a continuous modification of Z. We will choose 
such continuous version of Z from now on. 

If (O, J 7 , P), {J 7 *}, and (Y, Z, W) is one weak solution of (|1.7[) . we can con- 
sider the image measure of P under the mapping (Y, Z, W) : ft — > W 1+d x W d 
defined by 

u h-> ((Y(Q),Z(Q)),W(u)) , for we ft, 

which is denoted by Q. Since the projection of Q on the third component W is 
a Wiener measure on (W d , B w ), denoted by Q* from now on, and all the spaces 
are Polish under the uniform topology, there exists a unique regular conditional 
probability Q{-|w} such that: 

(1) for uj e B w , Q{-|w} is a probability measure on (W 1+d , B Y ® B z ); 

(2) for Ae B Y (g> B z ', the map lo ^ Q{A|cj} is S^-measurable; 

(3) for A e B Y ® B z and B e B w , we have 

Q(A x B) = [ Q{A\u}Q*(du;). 
Jb 

Definition 13 The weak solution to JJ.7| ) is called unique in law if for any two 
weak solutions (0, J",P), {F t } , (Y,Z,W) and (0, J",P), {T t }, {Y,Z,W), the 
probability distributions of (Y, Z) and (Y,Z) are equal, i.e. P(y,Z) = P(y z)- 

The weak solution to {1. 7| j is called pathwise unique if for any two weak 
solutions (Y, Z) and (Y,Z) defined on the same probability space (ft, J 7 , P) with 
the filtration {J 7 *} and the same Brownian motion W , (Y, Z) is a continuous 
modification of(Y,Z). 
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For a SDE, the celebrated Yamada-Watanabe theorem states that the weak 
existence and pathwise uniqueness of the solutions to a SDE implies the ex- 
istence of a strong solution. As Kurtz [53] pointed out: strong solution is a 
consequence of measurable selection, and such result has little to do with the 
equation, but really a consequence of the convexity of collections of the proba- 
bility distributions of solutions. If the compatibility constraint (Sec Definition 
^ is satisfied, we further have the adapteness of solutions. 

Theorem 14 If the coefficients satisfy Condition^ then the weak solution to 
BSDE J1.7] ) is pathwise unique, and the strong solution also exists. 

Wc first prove the pathwise uniqueness. By employing the Girsanov's theo- 
rem reversely, we have the following pathwise uniqueness result. 

Lemma 15 If the coefficients satisfy Condition fj[ then the weak solution to 
\1. 7| ) is pathwise unique. 

Proof. Suppose (Y, Z) and (Y, Z) are two weak solutions on (fi, J 7 , P) with 
{Ft] and Brownian motion W . By applying Ito's formula to e at (Y t — Y t ) 2 for 
some a to be determined, wc obtain 

e at {Y t -Y t ) 2 

= - 2 J e as (Y s - Y s )d(Y s -Y s )- J e as d[Y — Y,Y — Y] s 
-J ae as (Y s -Y s ) 2 ds 

£ e as (Y s - Y s ) ^J2(Zlf(s, ZD - Zif(s, Z{)) + (h(s, Z s ) - h(s, 2.)) j ds 

2 [' e as (Y s - Y.) YjZI - Z\)dWi - f e as £ \Z{ - Z\\ 2 ds 

ae as {Y s -Y s fds. (3.1) 



= 2 



Note that for s <G [0, T], and z j , z j € R for j = 1, • • • , d, 

\z j f j (s,z j ) - z j f j (s,z j )\ 

< \z J f(s, z j ) - z ] f j {s, z j )\ + \z j f j {s, z j ) - ?f (s, 

< d(T+ - z j \ + C^Wz 1 - z j \ 

< d{T+\z j \ + \zi\)\z j -z j \. 

Now if we set 

z s — z s 

when Z\ — Z J S ^ 0, and ft 8 = for s E [0, T] otherwise, then there exists a 
constant C 8 such that \ft s \ 2 < C 8 (l + \Z{\ 2 + \Z{\ 2 ). Based on such ft, we 
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define a new probability measure Q by fgj = £(N) where 

d „. 

N = zZ ti dW i> 

and under Q define a new Brownian motion B by B — W — [W, iV] . Then under 
the probability measure Q, (|3.1[) reduces to 

e Qt (r t -r t ) 2 

= -2 e as (Y s - Y s ) J2(Zi - /.: hUV + 2 / e as (Y s - Y s )(h(s, Z s ) - h(s, Z s ))ds 
Jt J=1 Jt 

j-T d t 

e as Y^\Zl- Zi\ 2 ds- ae as (Y s -Y s ) 2 ds. 
Jt j=1 Jt 

By taking expectation under Q we have 

E<Z[e at {Y t - Y t f] = E<* j e as 2(Y s - Y s )(h(s, Z s ) - h(s, Z s ))ds^ 

-E Q \J T e as \Z s - Z s \ 2 ds \- eQ \[ <*e a "(Y s " Y sfds \ ■ 



By the elementary inequality 2ab < A 2 a 2 + b 2 /X 2 

- VA 2 j. 

A 2 



2(Y S - n)(/i(s, Z s ) - h(s, Z s )) < X 2 (Y S - Ys) 2 + %\Z S - Z s \ 2 . 



By choosing A 2 = a and a = 2C 2 , we obtain 

so Y t = F t for t G [0,T], a.s., and Z t = Z t for a.e. i g [0,T], a.s.. Now the only 
step left is to verify S'(N) is a uniformly-integrablc martingale, and we need to 
verify N is a i?AfO-martingale under P. In fact for any stopping time t <T, 

E{\N T - N T \ 2 \T T } = E ^ l^l 2 ^} 

<C 8 E^ (d+\Z s \ 2 + \Z s \ 2 )d,s\T T ^, a.s.. 

The way to control the integral term involving Z and Z has already been pre- 
sented in the proof of Lemma [T2] Therefore Q defined above is indeed a prob- 
ability measure. ■ 

The idea for the following lemma is standard: to transfer the structure of a 
weak solution such that W 1+d becomes the sample path space for (Y, Z) and 
W d that for W. What allows us to carry it through is the regular conditional 
probability introduced above. 
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Lemma 16 For BSDE |i.7| j with the coefficients satisfying Condition^ path- 
wise uniqueness implies uniqueness in law, and furthermore if the weak solution 
exists, then the strong solution also exists. 



Proof. Since the proof is quite standard, we only present the basic steps. 

Let (n,F,P), {F t }, (Y,Z,W) and (Cl,F,P), {F t }, (Y,Z,W) be two weak 
solutions. Let W 1+d and W 1+d be two copies of C([0,T]jii) x C([0, T];R d ). 
By using the regular conditional probability Q{-|w} and Q{-|w}, we define a 
probability measure tt on the probability space (9,6(6)) by 

tt {(dy, dz), (dy, dz),diu) = Q{(dy, dz)\uj}Q{(dy, dz)\uj}Q* (du) 

where 

(6, B(Q)) = (W 1+d x W 1+d x W d , B Y ® B z ® B Y <g> ® B w ). 

On (9, 23(9), 7r), we further define the filtration {Q t } which is generated by 
c((y s , z s ), (y s , z s ),lu s : s < t) augmented by the 7r-null sets in £>(9). Then under 
tt and {Qt}, w is still a Brownian motion. In fact by the compatibility constraint 
in Definition 3, B Y <g> Bf is independent of Bf . Hence for A t £ Bj <g> Bf , 

Q{A t |w} = Q{^|w t }. 

Likewise we also have Q{A t |w} = Q{A t |w t } for A t £ Bj ®Bf . Based on above 
relationships, for B t £ Bf , u £ [t,T] and £ € 



Q{A t | w }Q{A t | w }Q*(da 



B 4 



D i(£,Wu-a>t)d/-|*/ 



Q*(dw) / Q{A t |a; t }Q{A t |w t }Q*(dw) 



B, 



E* [IaJaM] 



which means {lo u — Lo t : t < u < T} is independent of the filtration {Qt}- 

Therefore (y, z,uS) and (y, z,ui) are two weak solutions on the same filtered 
probability space (9, B(<3), {Qt}, i"). Pathwise uniqueness means 



* ({(y, z )i (v, z )> u} € e : (y, z) = (y, z)) = 1, 



(3.2) 



so for any A £ i3 r <g)£> z , the probability distribution P ( y Z )(u; e 17 : (Y, Z) e A) 
equals 

7r({(y,z),(y,z),w} £ 9 : (y, z) € A) = tt ({(y, z), (y, z), w} £ 9 : (y, z) e A) 

which is equal to the probability distribution P(y- ^(cj <G Cl : (Y, Z) £ A). 

To prove the second claim, we firstly show Q{-|oj} and Q{-|w} assign full 
measure to the same singleton. By (|3.2p and the definition of tt, we have 



W" J(y,z) = (y,z) 



Q{(dy,dz)\u}Q{(y,z)\u}Q*(du) = 1, 
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so there exists N G B w with Q*(N) = such that 

f Q{(dy,dz)\u}Q{(y,z)\uj} = l, for^r. 

This can only occur if there exists a S^/2? y <g>£> z -measurable map $ = ($ y , $ z ) : 
W d -> W 1+d such that 

Q{{y,z)\u} = Q{(y,z)\u} = 5 9{u) (y,z), for u G N c . 

It then follows (y, z) = ($ y (w), $ z (w)) for we]V c . But recalling w n- Q{A t |w} 
is S^-measurablc for At G (g) £>f , then by the standard Dynkin arguments, 
$ is in fact also BY /Bj ® £> t z -mcasurable, and on any given filtered probability 
space (fi, J", J" t ,P) satisfying the usual conditions with W being a Brownian 
motion on it, 

$ y (W t ) = 0(W T ) + ^ ft(a, $ z (Wg)ds 

+ [ Y,® Z ' j {W s )P{s,® z > j {W s ))ds- [ J2^ j (W s )dWi, 
so $ is a strong solution. ■ 



4 Optimal portfolio in incomplete markets 

In this section our main aim is to demonstrate how the functional differential 
equation approach and the nonlinear Girsanov's transformation can be used in 
finance. Specifically we consider an example of optimal portfolio problems in 
incomplete markets which is often used in the indifference valuation. 

Assumption 1 (the probability space) Let (Q, J- , P) be a complete probability 
space which is to be determined, and {J^t} be its associated filtration satisfying 
the usual conditions, which is also to be determined. 

Assumption 2 (the market) The market is built with three assets: a risk- free 
bond with zero interest rate, a tradeable asset and a nontradeable asset. The 
pricing dynamic of the tradeable asset satisfies the following SDE on the above 
given probability space: 

r ds t /s t = fifdt + afdWt, (41) 

\ So — s, 

and the dynamic of the nontradeable asset follows 

f dV t /V t =vYdt + aYdW t v + aYdW t , . . 

\V = v 

on the same given probability space, where W = (W v , W) is a two-dimensional 
Brownian motion to be determined. a v and y, % ,a % for i = S,V, as the market 
coefficients, are bounded and continuous functions. 
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Assumption 3 (the investor) The investor has an exponential utility function 
depending on his/her terminal wealth, which has the form: 

U(x) = -e -7 *, for x e R, 

where 7 > representing the degree of the investor's risk aversion. 

Assumption 4 (the trading strategy) The investor, with initial wealth x, in- 
vests in the tradeable asset and the risk-free bond during the time period [0, T] . 
Let 7r be the amount of money invested in the tradeable asset. We assume tt is 
taken from the following admissible set, which of course depends on the above 
given probability space. 

■Aad '■= : [0,T]xfi — > R : tt is J- \- adapted, self -financing and |lvr|[jy2 ro,T] < 
The dynamic of the investor's wealth process, denoted by X x {tt), follows 

dX?(ir)=rr t (rfdt + afdW t ), 
X*{tt) = x. 

Assumption 5 (the cost functional) At time t = T, the investor gets the total 
amount X^(tt) plus a random endowment g(Vr, St), where g is assumed to be 
Lipschitz continuous and uniform bounded. The investor decides the optimal 
trading strategy to maximize the following cost functional: 

Here we use the superscript P to emphasize the expectation is taken under the 
probability measure P, which is to be determined. 

The random endowment g depends not only on the nontradeable asset V 
but also on the tradeable asset S, which distinguishes the current problem from 
the ones usually considered in the literature. In the book edited by Carmona 
[PJ, this problem is even called an open problem (|19j). Such form of random 
endowment actually appears naturally when one wants to consider the credit 
risk of options traded in OTC markets (see Henderson and Liang [30] ). We also 
emphasize the well known Cole-Hopf transformation does not help to deduce 
the closed-from solutions in our setting. 

In the following we give the definition of weak admissible trading strategy 
and the corresponding weak formulation of optimal portfolio problems. For the 
weak formulation of general stochastic control problems, we refer to Yong and 
Zhou |H]. 

Definition 17 A triple (0, P) {J-t} and (tt, W) is called a weak admissible 
trading strategy if 

(1) (f2, J 7 , P) is a complete probability space with the filtration {J 7 *} satisfying 
the usual conditions; 

(2) W is a Brownian motion, and the increment {W u — W f : t <u < T} must 
be independent of a -algebra Tt; 

(3) tt is taken from the admissible set A ac i- 



sup E p 
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The set of all weak admissible trading strategies is denoted as A^ d , and a 
generic element in such weak admissible set A^ d is denoted as II. The investor 
decides the optimal weak admissible trading strategy II in order to maximize 
his/her cost junctional: 



sup E 



_ e -~t(X*(TT)+g(VT,S T )) 



(4.4) 



Remark 18 The motivation of introducing the above weak formulation of op- 
timal portfolio problems is more from mathematics rather than finance. Later 
We will use the martingale optimality principle to deduce a quadratic BSDE 
as the characterization of the optimal portfolio, and we will look for the weak 
solution of such quadratic BSDE. The probability space will be chosen from the 
weak solution of the associated quadratic BSDE. 

Next we use the martingale optimality principle to characterize the opti- 
mal portfolio. For a given filtered probability space (fi, F, Ft, P) with a two- 
dimensional Brownian motion W = (W v , W), all of which are to be determined, 
we want to construct a family of stochastic processes 



( 



-y(X?(iv)+Y t ) 



te[o,T] 



such that 
(1) the process (- 



- 7 p<: t *M+y t ) 



)te[o,T] is a supermartingale for any 7r G A a d, 
and there exists an optimal n* £ A a d such that (— e _7 ( Jf t'( 7r ^ +Yt ')te[o,T] is a 
martingale; 

(2) the auxiliary process (^)te[o,T] has the terminal value Yt = g(VT,Sr)- 
If such auxiliary process (y*)te[o,T] an d the optimal 7r* exist, then we have 

E P \- e -~f( X TW+Y T )] < _ e -7(*+>o) ; for any n g Aadi 



and 

E 1 

Therefore 



- 7 (XJ(7r*)+r T ) 



= - e -^ x+Yo \ for optimal tt* e A ad . 



sup E 



3 - 7 (x?( OT )+y T ) 



E* 



--f(X*(n*)+YT) 



-y(x+Y ) 



Note that the filtered probability space (f2, F, Ft, P) and the Brownian mo- 
tion W are still to be determined. Next we use the weak solution of a quadratic 
BSDE to characterize the auxiliary processes Y and 7r* , which also provides us 
with the filtered probability space (il, F, Ft, P) and the Brownian motion W. 

Theorem 19 Let (0,J",P), {F t } and (Y,Z,W) (with Z = (Z V ,Z)) be the 
weak solution to the following quadratic BSDE: 

Y t =g(V T ,S T )- J feds- J (Z v dWY + Z s dW s ) (4.5) 
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with 



2^' 2 7 (af)2' 
TTien the value function of the optimal portfolio problem * s given by 

and the optimal weak admissible trading strategy IT is the triple (0, J- , P), {Ft} 
and (77* , W) with 

<=-Zt + -{^. (4.6) 
IV? t ) 2 

Proof. On a given filtered probability space F, Ft, P) with the Brownian 
motion W, we suppose iYt) t ^[o,T] satisfies the following BSDE: 

Y t = g{V T ,S T )- J fsds- J (zYdWY + Z s dW s ), 

where the driver / is to be determined. By applying Ito's formula to e~ J ^ x t^ +Yt \ 
we obtain 



Vl 1 ! -I- 1 
7 2 



2 



= e->W«+*> |-7^ - 7/* + y [(^f) 2 - t 2 + (^) 2 + (^) 2 
+2(rf Z t 7r t ] } + martingale term. 

Since (— e _7 ^ J,<: * : ' 7r ^ +Y ^) te [o.T] is supermartingale for any n € _4 a( 2, and a mar- 
tingale for optimal 7r* € -4 a£ i, we must have and 7r* such that 

— (fi t ) [n t + Z t +— / t + -rg-Zt - s - "I ft > 

for any 7r G -A a d, and equality holds for optimal tt*. Therefore by solving the 
above variational inequality, we obtain 

t — T ' (vV\2 i Mt ^ (Mt ) 2 
2 of 27(crf) 2 

and 

7T t = -Zi 



7(of) 2 ' 
■ 

Iii the following we will employ the functional differential equation approach 
and the nonlinear Girsanov's transformation to find the weak solution of BSDE 
(|4.5[) . Since the coefficients satisfy Condition [3J by Theorem [Til the weak so- 
lution we will find is pathwise unique, and moreover the strong solution also 
exists. 
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The idea is to use the strong solution of the following FBSDE (|4.7|) to con- 
struct the weak solution of BSDE (|4.5|) . Let's start with a Brownian motion 
B = (B V ,B) on (ft, J 7 , Q) with the filtration {J 7 *} satisfying the usual condi- 
tions, and consider the following FBSDE: 



dkiVt 



1 



-aY^dt + aYdBY 



07^, 



InFn = 



din 
In So 

dY t 
Yr 



lav 
1 



-((t?ydt + it?dB u 



(4.7) 



= In ; 



2 7 (af) 2 



+ Z t v dB t y + Z t dB t 



g(e laVT ,e lnST ). 



Note that FBSDE (j4.7[) is linear and the coefficients satisfy Condition [TJ so by 
Lemma[5J we know there exists a unique solution (Y, Z, In V, In 5) € C([0, T]; i?) x 
if 2 ( [0 , T] ; i? 2 ) x C ( [0, T] ; JJ) x C ( [0, T] ; J2) . 

Based on the solution (Y, Z), we define a new probability measure P by 



dP 
dQ 



S{N), 



where $ (N) is the Doleans-Dade exponential of N with 



N 



IzYdB] 



^dB t . 
of 



By Lemma fTJl we know P is indeed a probability measure. Under the new 
probability measure P, by Girsanov's theorem, W = B — [B, N] is a Brownian 
motion with 

' W V = B V 



-ZYdt, 

z 



W = B - 



%dt. 

o of 



Under the probability measure P and with the Brownian motion W, let's 
rewrite the backward equation in FBSDE (|4.7|l : 



dYt 



(/40 



S\2 



27(of) 



dW i 



v , 7 ^y. 



~z; dt) +Z t dW t + ^dt 



i4 



= f t dt + zYdw t v +z t dw t 

with Yt = g(Vr, St), and rewrite the forward equations in FBSDE (|4.7[) : 



dVt/Vt = \i\dt + oYdWY + a(dWt 
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with 

dSt/S t = t^dt + afdWt. 

Therefore the triple (Q, J- , P), {J~t} and (Y, Z, W) is just one weak solution we 
want to find. 
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